A control protocol of finite dimensional quantum systems 
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An exact and analytic control protocol of two types of finite dimensional quantum systems is 
proposed. The system can be drive to an arbitrary target state using cosine classical fields in finite 
cycles. The control parameters which are time periods of interaction between systems and control 
fields in each cycles are connected with the probability amplitudes of target states via triangular 
functions and can be determined analytically. 
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Quantum control is to drive a quantum system from 
an initial state to an arbitrary target state through its in- 
teraction with classical control fields or with a quantum 
accessor. It was first proposed by Huang et. al. in 1983 1 1] 
and then attracted much attention of chemists, physicists 
and control scientists. Various notations in classical con- 
trol theory were generalized to the quantum control, such 
as open and closed control, optimal control [2| , controlla- 
bility feedback control [3] and so on. Coherent and 
incoherent (indirect) control schemes are proposed. In 
later case the system is controlled by its interaction with 
a quantum accessor which is controlled by classical fields 
[8H11| . Typically, in the approach of quantum control, 
one should first model the controlled system and exam- 
ine its controllability which is determined by the system 
Hamiltonian and interaction Hamiltonian with classical 
fields, and then design classical fields to stream the sys- 
tem to the given target state, which is referred to as 
the control protocol and is the issue we would like to ad- 
dress in this paper. Some works were proposed along this 
line, for example, using the Cartan decomposition of Lie 
groups 

In this paper we shall develop an explicit control proto- 
col of finite quantum system with (1) all distinct energy 
gaps such the Mores potential, and (2) all equal energy 
gaps except the first one. We use cosine classical field 
to drive the quantum system to arbitrary target states 
in finite cycles and the control parameters are interac- 
tion time intervals between system and control field in 
each cycles. Control parameters are linked with proba- 
bility amplitudes of the target states through triangular 
functions and can be obtained analytically. 

This paper is organized as follows. In Sec. II, we for- 
mulate the controlled system and control scheme and in- 
vestigated the controllability. In Sec. Ill, we present the 
control protocol of system with all distinct energy gaps 
and in Sec. IV we consider the system with equal energy 
gaps except one. We conclude in Sec.V. 



A. Control Systems 

Consider an V-dimensional non-degenerate quantum 
system with eigen energy E n and corresponding eigen- 
states |n), described by the Hamiltonian 



iV 



H = J2E n \n)(n\ 



(1) 



n=l 



Without losing generality, we assume Hq is traceless, 
namely ItHq = 0. In this paper, we only consider two 
different types of systems, the first one having all equal 
energy gaps except the first one, namely 



Mi A*2 = Ma = • • • = Miv-i) 



(2) 



where fa = Ei+x — Ei is the energy gap, and another one 
with all district energy gaps 



Hi^fMj, i^j = 1,2, 



,N- 1, 



(3) 



For later convenience, we call them the system I and sys- 
tem II, respectively. For System I, we also define energy 
gaps 

fa = E i+1 -E 1 , i = l,2,.--,JV-l. 

The purpose of this paper is to develop a control 
scheme to drive the systems to an arbitrary target state 
from an initial state, using some independent classical 
fields f m (t)- The total Hamiltonian of the system and 
control fields can be generally written as 



M 



H 



Hq + Hi, Hi — f m (t)H n 



(4) 
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where M is the number of independent classical fields. 

For the TV-dimensional systems considered in this pa- 
per, the total control process includes N — 1 cycles. In 
the m-th cycle, we first apply a classical field 

fm(t) = £ m cos(v m t), 

where v m is the frequency, to control the system for 
time period r m , and then turn off the control field such 
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FIG. 1. Control fields. 
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that the system evaluates freely for a time period T' m , as 
showed in FigfT] 

For system I, the frequency of the control field is chosen 

as 



&m = (Em+1 — Ei) k 



(5) 



It is easy to find that u m ^ uj n for m ^ n. This means in 
each cycle, only transition between E m and E\ occurs. 
For system II, the frequency v m is chosen as 



k>m — {E m +i — E m ) h 



(6) 



This means that, in each cycle, only transition between 
level E m and E m+ i occurs, as all u m are different. So 
the control process includes N — 1 cycles. 

For both systems, there are two processes in each cycle. 
We first apply the control field to control the system for 
a time period r m , and then turn off the control field and 
allow the system to evaluate for a time period T' m . We 
will see that the first process provides the real probabil- 
ity amplitude of the target state and the second process 
provides the phases. Therefore, the control field can be 
rewritten as 



fm(t) = 



£ m cos(V m i), t m -i < t < t m -i 
0, otherwise, 



(7) 



where t m — X)fcLi( T /c+ T jc)- Those N— 1 control fields are 
independent in the sense that each f m (t) ^ in different 
time period. 

The whole control process can be cquivalcntly regarded 
as control by one control field f(t) = J2m=i fmif), where 
fit) is shown in Fig.l. 

For system I, in the m-th cycle, we apply the classical 
field 



f {m \t) =U cos (v m t), v„ 



(8) 



which causes transition between energy levels E\ and 
E m+ \. So the interaction Hamiltonian between control 
field and system in Schroding picture is 

H im) =/(«•) ( t ) flm (|l)(m + 1| + \m+ 1)<1|) 

=n m cos{v m t) (|l)(m + 1| + \m + 1)(1|) , (9) 

where 

{Jm is the coupling constant and — £,m9m- 
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FIG. 2. Interaction between systems and control fields in each 
cycle. For system I, control field causes transition between the 
level Ei and level E m +i, while for system II, the control field 
causes transition between adjacent energy levels. 



For system II, the interaction Hamiltonian between the 
system and the control field is 

H (m) =/ (m) (t)ffm ( | m) ( m + 1| + | m + 1) ( m |) 

=ft m cos(V m t) (|to) (m+ 1| + h.c) , (10) 
where g m and Q m = Cm9m as in system I, but frequency 

Interaction between systems and controls are illus- 
trated in Fig.©. 



B. Complete controllability 

Before presentation of the control protocol, we first ex- 
amine the controllability of this control scheme, namely, 
to examine whether the Lie algebra generated by the 
skew-Hermitian operators iHq and iH m 



C = Gen{iH ,iH m \m = 1, 2, N - 1} 

is su(iV) Q. Here H m for system I and II are 

H m = |l)(m + l| + |m + l)(l|, 
H m = \m)(m + 1| + \m + 



(11) 



(12) 
(13) 



respectively. Or equivalently, iHo,iH m generate the 
Chevalley basis of su(A^) |, [| 



ix n — i(\n)(n + 1| + \n + l)(n|), 
iy n = \n)(n + 1| - |n+ l)(n| , 
ih n — i(\n)(n\ — \n + l){n + 1|), 



(14) 



where n = 1, 2, N — 1. In fact, it is enough to prove 
ix n £ C (or iy n £ £), as iy n = ^[iHa^Xn] and ih n = 
-[ix n ,iy n ]/2. 



3 



For system I, it is obvious that 

ix 1 =iH 1 = i(\l)(2\ + \2)(l\)eC. (15) 

Then 

iy 2 = [iiT 2 ,iHi] = |2)(3| - )3> <2| e £, 

ix 2 = —[[iH2,m],iH Q ] = i(|2){3| + |3)(2|) e C. (16) 



Recursively, we have 

iy™, = Mm* [[iH m ,iH m -x],iHo] ■ 

So the system I is completely controllable. 
For System II, iH m itself 



(17) 



iH m = i(\m)(m + 1| + |m+ l)(m|), (18) 

is nothing but the generator ia; m . Therefore, the system 
II is completely controllable. 



III. CONTROL PROTOCOL OF SYSTEM I 

In this section we investigate the control protocol of 
System I. Suppose the system is initially on the ground 
state |^o ) = |1)- The system is driven to an arbitrary 
target state after N—l cycles. Relationship between the 
control parameters {T m ,r^j} and probability amplitude 
of target states is explicitly established. 



A. Interaction Hamiltonian 

Changing to the interaction picture, we obtain 

H { j m) =C7 f (t) [n m coa(p m t) 

x (|l)(m + l| + |m + l)(l|)]Z7 (t) 

a, 



J(v m — U m )t 



1} (m + 1 



+ e -*("m+Om)t |^ ^ + J_| 

+e^-^)*|m + l)(l|) , 



(19) 



where Uo(t) = e - lH ot/h^ anc j we nave usec j cos(i/ m t) = 
( e w m t _)_ e ^ lu rnty2. As we require that the control field 
is resonant with the levels E\ and E m , namely, v m = 
uj m , we can neglect the high-oscillating terms e ±i ^ m+Vm ' t 
under the rotating wave approximation [131 ] . We finally 
obtain 



tf} m) = ^(|l><m + l| + |m + l>(l|), (20) 

which does not depend on time t explicitly. So the time 
evolution operator in interaction picture can be written 

as 



u*r\t) = e- iH( r h i n . 



From the fact 
n > 0, 



^ 2n 






2 



(|l)(l| + |m + l)(m + l|) 



^ 2n+l 


i2 m 




2 



2n+l 



(|l)(m+l|+h.c), 

J 

n > 0, 

we obtain the time evolution operator as 

U\ m \t) =1+ (cos(n^t) ~ 1) (|1)(1| + \m + l)(m + 1|) 
-tsin(fi^t)(|l)(m+l|+h.c), (22) 

where Q,' m — Q m /2h. 

B. Cycle 1 

Suppose that the system is initially on the state |1). 
When interacting with control field for time period n, 
the system is on the state 

1^=^)11), 

= cos(nin) |1) - isin(niri) |2) . (23) 

Changing back to the Schroding picture, we have 

\y Jl ) s =e-* E ^/ h cos(n' 1 T 1 )\l) 



le 



-iEiTt/h 



sin(fiin) |2) . (24) 



We then turn off the external field and allow the system 
to evaluate for time period t[. We get 



^ = e -^ rl/H = (1) | 1} + fl (D | 2) ^ (25) 



with 



cos(r2' 1 Ti), 



4 1} = - e - tE2iTl+ < )h isin(fiin). 
C. Cycle 2 



(26) 



For the cycle 2, the initial state is the the final state of 
cycle 1, namely, the state (1251) . We first apply the control 
field f2{t) — £ 2 cos(^ 2 i) for time period r 2 . Using (|2"2"|) for 
to = 2, we obtain the state in the interaction picture 

=^i(r 2 ) |Vi>'/ 

= cos(0 2 t 2 )4 1) |1) + 4 X) |2) - i sin(n 2 T 2 )a[ 1) |3) . 

(27) 

Changing back to Schrodinger picture, and after free evo- 
lution for time period r 2 , the state is 



(21) 



\<h)s =C/ o(r 2 ) |V 2 ) S = (7o(t 2 )(7 (t 2 ) \i) 2 )j 

=ap> [1) + 4 2) |2> + 4 2) |3) , (28) 



4 



where 



(2) _ lEl (T 2 +r: 2 )/h 



cos(f2' 2 T2) a i 



(1) 



a (2) =e - iE2 (r 2 +^)/ fia (l) ) 



4 2) = -e- lE3 ^+^^ h tsm(n' 2 T 2 )a ( ^ . (29) 



D. From (m — l)-th to m-th cycle 

To obtain the explicit expression of the target state, 
we need the recursion relation between coefficients of the 
m-th cycle and (m — l)-th cycle. Suppose that after the 
(rn — l)-th cycle, we obtain the state 



\ip m -i)' s = 514™ 1} \ k ) 



(30) 



k=l 



Interacting with control field for time period r m , we have 
the state in Schroding picture 

H m ) s = e- lE ^/ h cos(n' m T m )a { ™- 1} |1) 

m 

+ J2 a i m ~ 1)e ~ lEkTm/h \ k ) 

2 

iEm +1 r m /n isin ^ mTm)a (rn-l) ^ + ^ _ 



fc=2 

e~ 



After free evolution for time period r' m , the final state of 
the m-th cycle is 



m+l 



\iPm)' s = e~ iH ^ n \^ m ) s = 4 m) l fc ) . (3 1 ) 



fc=l 



where the coefficients are 



M =p -^T m /, cos( ^ Tm)a ( m -i 



a; = e 



(32) 



4 m) = e - iE » T ™l n ai m - 1} , 2 < k < to, (33) 



(m) 
Z m+1 



^ +1 ^/n isin(n ^ Tm)a (m-i) 



and we use the notation 



(34) 



(35) 



hereafter. Eqs. (|32H3~4"]) establishes the relationship be- 
tween the probability amplitudes of the (m — l)-th cycle 
and the m-th cycle. 



E. Target state 

For 2 and 3 dimensional system, the target state has 
been found in subsection B and C. So we suppose that 
N > 4 in the rest of this section. From (|3"2"|) and ([2"fJf , we 
can easily find that 



(m) 

a\ = exp 



iEi 



h ^ 

i=l 



[Jcosffi^). (36) 



As a™li depends on a^ m L> only, we can easily find that 



(m-l) 



(m) 
1 m+\ 



6XP \~h 



E m +iT m + Ei 53 Ti 

7=1 

771—1 

sin(fi^r m ) Y[ cos(n-Ti). 



(37) 



For other coefficients, using (|29p . we obtain the explicit 
probability amplitude 



,( m ) 



,( m ) 



exp 



exp ■ 



iE? 



5> 



i=i 



i sin(J7' 1 Ti), 



E 3 J2 T * + E ^ 



(38) 



To derive coefficient fflj. , 3 < /c < to, we use 



7 (m-l) 



6XP 



X i 



E m T m -i + Ei 53 ^ 

i=l 

7n-2 

sin(fi^_ 1 r m _i) J] cosCfiJri), (39) 



which is obtained from (|37[) by replacing to by to — I. 
Then we can recursively have 



a (rn) =e -iWH a ^ =exp 



iE„ 



exp ■ 



i= m 

m m— 2 



-,m — 1 



7— 7n — 1 

777 -2 



7 = 1 



x isin^^.jTm-i) ]J cos(0/jTi). 

7=1 

According to (f3T)|) and (|4"0|) . we can obtain a£ 

m _ p -iE k T m /H n {m-l) 



— e -iE h T m /n e -iB k T m -i/h a (m-2) 



(40) 



(T m +T m _i+-+T* + i)-(fc) 



exp 



fe-2 



E k J] T + ^^T 



7 = fe-l 

fc-2 



X I 



ism(Q' k _iT k -i) Y[ cos(fi-Ti), 

7 = 1 

(3 < k < m), 



(41) 



where we have used the result of l ' given by ([3^ 
with m replaced by fe. Notice that when k = m + l, BT] 
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recovers the a!^ 1 given in ([571) . Therefore (HTj) is also 
valid for fc = m + 1, and all the probability amplitudes 
after m-th cycle are given by (j36]) . (|38|) and (|4Tj) with 
3 < fc< m + 1. 

For the system we considered here with dimension N, 
we need N — 1 cycles to arrive at arbitrary target states. 
Letting m = N — 1, we obtain the probability amplitude 
of the target state 



(N-l) 



(JV-1) 



: exp 



exp 



exp ■ 



N-l ' 

i=l 
7V-1 



Af-1 



n cos (^^) ) 
i 

i sin(ri' 1 ri), 



N-l 



fc-2 



=k-i 

fe-2 



x isin^.iTfe-i) ]J cos(fi-Ti), 
»=i 

(3 < fc < AT). (42) 



For a given target state, namely, C n and 7„ are 
given, we can calculate control parameters {r„,r^|n = 
1, 2, N — 1}. From (l45j) and C2, we can determine n. 
Then form (1461) with n = 3, we can obtain T2 from C3. 
Repeating this process, we can obtain all parameters r n , 
ri = 1,2,..., N- 1 from (gSJ). 

All t[ can be obtained from (|49p. From 72 and 73, we 
can obtain 



JV-l 



E r - 



AT-1 

1=2 



from which we find T\ and X^il? 1 ^i- From 74, we find 

JV-l 

S 4 ^T J + ^ 1 (T 1 +T 2 ), (51) 



EiT lt 



(50) 



i=3 



from which as well as T\ and yi— o Tj, we can obtain T 2 . 
Repeating this process, we can obtain all Ti and thus all 



F. Control parameters 

For a control problem, the target state, or in other 
words, the amplitude a^ -1 of the target state, is given. 
What we need to do is to determine the control param- 
eters { , I i — 1,2, ...,N — 1} from the probability 
amplitude of the target state. For convenience, we write 
the target state as 



JV 



W) = ^2lnC n \n) 



(43) 



71=1 

where C n s are the real part of the amplitude 

d = J] cosfe), (44) 

i=l 

C 2 = sin(niTi), (45) 
n-2 

C n = siii^jrn-i) [| cos(fijTi) , 3 < n < N, (46) 

i=l 

and 7„s are phases 



71 = exp 



72 



7n 



i exp 



JV-l 

»=i 

JV-l 

~ir5> 



i = l 



= -zexp|-- 
(3 < n < N). 



N-l 



n-2 



(47) 
(48) 

(49) 



IV. CONTROL PROTOCOL OF SYSTEM II 

A. Time evolution operator 

For system II, the interaction Hamiltonian is given in 
Ea. (fT0)) . This Hamiltonian is same as ([9]) for system I 
except the state |1) is replaced by \m). So we can follow 
exactly the same procedure as in last section, namely, 
changing to the interaction picture, using rotating wave 
approximation, and obtaining a time-independent Hamil- 
tonian in interaction picture 



H ,{m) = On ()m) ( m + 1| + | m + 1) ( m |) . (52) 



Using 



(H>(m)y n = (| m ) ( m | + \ m + l) (m + l\) 



(n > 0) 



2n+l / Q 



2n+l 



(|m) (m+l\ 



1 1 2 

+ \m + l)(m\), (n>0), 

we can obtain the time evolution operator in the inter- 
action picture 

C/j m V) = / + [cos(O^f) - l](\m)(m\ + \m + l)(m + 1|) 
-ism(il' m t)(\m){m + l\ + \m + l)(m\), (53) 

where Q' m = Q m /(2h). 
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B. Determine amplitude a m 

For this model, the cycle 1 is exactly the same as the 
system I. So after the cycle 1, the system is driven to the 
state 



where 



l^i>s = oi 1) H>+4 1) |2>. 



i« = e - iE ^/ % cos(fiiri), 



i sin(S7^ ti). 



(54) 



(55) 



Different from system I, in cycle 2, the control field 
f(t) = £2 cos(w2i) causes transition between |2) and |3). 
We can find the state after the cycle 2 in Schrodinger 
picture as 



|^ 2 )' s = 4 2) |f) + 4 2) |2)+a^|3), 



,(2) 



where 



fl (2) =e -iBiT a //l a (D j 
a (2) =e -iE 2 T 2 / ftcos(0 / T2)a (l) 

4 2) =-e-^ r2 /^sin(0^r 2 )4 1) . 

To obtain the target state, we first find the recursion 
relations between the (to— l)-th cycle and the m-th cycle. 
To this end, we suppose that, after m — 1 cycles, the 
system is on the state 



m-l/s 



E< 

fe=i 



(m-l) 



(60) 



Then after interactions with the control field for time 
period r m , and free evolution for time period r' m , we find 
the final state after cycle to as 



m+i 



k=l 

with (m > 2) 

a M = e -iB k r m /ft a (Ti»-i) ) 1 < fc < m _ i, 
a M = e -«V»r m /» cos( jy m7Vn)o O r -l) > 

ism(^r m )a^- 1 ). 



_ _ p - 1 £ ra+ iT ra /S 



(61) 



(62) 
(63) 
(64) 



"?7l + l 

From those recursion relations, and initial conditions (|55[) 
we can find all the explicit expressions of . It is easy 
to see that 

J™) 



exp 



exp 



i=l 
■ zt 1 m 

ft ^ 

t=i 



cos(f2iri), (65) 
T; cos(r2^r 2 )isin(rj / 1 Ti), (66) 



(m) 



-iE m+1 T m /h 



/ .\m—l 

(—1) exp 



isin(f2^T m ) 



,(m-l) 



(1) 
2 



: exp 



lit 

- % -^E i+l Ti-i^ 



£J sin(ft -r^a 

t=2 
n 

nsin(njTi). (67) 



Then using (1521 and a m ™ 1 , which is obtained from (|6"7} 
by replacing to by to — 1, we have 

„M — p -iE k (T m W m )/h Jrn-l) 
u k ~ c a k 

_ p -iE k { Tm W m )/n -iE k {r m - 1 W m _ 1 )/n Jjn-2) 



(56) 


= exp 


iE k 

h 


m 

E T * 




= exp 


h 


m 

E T * 


(57) 








(58) 


= exp 


iE k 

h 


m 


(59) 






i—k 



iE k T k h 



cos(fl' k T k )a k ' 



(fc-i) 



k-l 



= 1 

fc-1 



*f (k - 1) cos(f7' fe r /c ) sin(f^Ti 



(68) 



i=l 



One can check that ([55)) includes the case fc = 2 and 
k = to as special cases. 

Therefore, after AT — 1 cycles, we arrive at the target 
state 



(N-l) 



a 



(N-l) 



■ exp 



: exp 



iE x 
iE„ 



N-l 



COS(OiTi) 



E^ 

i=l 

N-l . m-l 

m— 1 



^(to-1) cos(0^r m ) sin(Q-Ti), 



(2 < m < N- 1) 



N 



■ exp 



AT-l 



JV-l 



n sin fe). 



(69) 



C. Control parameters 

To determine control parameters rj, r/, 1 < i < — 1, 
we write the target state as 

iv>) = E a « i n ) = E 7nC " ' ( 7 °) 
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in which 



C\ = COS^iTi) 



C m = cos{Cl' m r m ) Y[ sin^rj), (2 < m < N - 1), 



i=l 



AT-1 



i=l 



(71) 



From 74, we can obtain 



e 4 ]T 



E3T2 



(74) 



i=3 



from which we obtain T2 as well as ^2^^ Ti. Repeating 
this process, we can obtain all and thus t[. 



CONCLUSION 



and phase 7„ 



71 = exp 



7 m = exp 



7w = exp 



' .„ N-l 



i=l 
7V-1 



. JV-1 
'ft 



^Si+xTi-i^iV-l) 



i=l 



(72) 



For a given target state, namely, C„ and 7„ are given, 
we can determine the control parameters {r„,T^|n = 
1, 2, iV— 1}. From C\ we can determine n, and then T2 
from C2 and the obtained t\ . Recursively we can obtain 
all r n . 

For r' n , from 72 and 73, we obtain 



JV-l 



N-l 



E ^" ^E + M ' 



(73) 



i=2 



respectively. As E 2 7^ S3, we obtain Ti and ^Z^ 1 -^- 



In this paper we proposed a protocol to drive two types 
of finite dimensional quantum system to an arbitrary 
given target states. The control parameters are time peri- 
ods {r m , T m \m — 1, 2, • ■ • ,N — 1} which can be explicitly 
determined from the probability amplitudes of the given 
target states. Relationship between control parameters 
and amplitudes is triangular functions and can be solved 
explicitly. The control fields in this protocol is the usual 
electric field described by the cosine function. 

We have 2(N— 1) real control parameters. In the target 
state there are TV complex or 2N real parameters. Taking 
into account the normalization condition of target state, 
one has 2(N —I) real parameters, the same as the number 
of the real control parameters. From this fact we can 
conclude that we can drive the system to an arbitrary 
target state by choosing appropriate control parameters 

{ T m, T m }- 

As further works, we would like to consider the indirect 
control protocol of finite quantum system by generalizing 
the control scheme in this paper. We also would like to 
consider the control protocol in the presentence of envi- 
ronment. 
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